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INTRODUCTION 


The  purpose  of  this  paper,  is  to  provide  a  means  of  estimating  and  comparing  the  wavelength 
conversion  efficiency  (figure  of  merit,  FOM)  of  new  nonlinear  materials  knowing  only  their 
bandgap.  The  approach  for  the  estimates  is  trend  analysis  of  existing  data '  for  the  second  order 
nonlinear  coefficient,  X®  or  dy  =  X^^^ijk/2 ,  the  index  of  refraction  n,  and  the  bandgap  E.  No 
underlying  physical  basis  for  the  result  will  be  presented.  The  result  of  this  trend  analysis,  in  the 
form  of  plots  of  x®  vs  E  and  FOM  vs  E  and  the  associated  fit  equations,  will  be  used  to  estimate 
X®  and  the  FOM  for  wavelength  conversion  efficiency  for  some  new  materials.  The  accuracy  goal 
agrees  within  a  factor  of  2-5  with  the  experimental  values,  which  is  quite  good  considering  the 
values  range  over  four  to  seven  orders  of  magnitude. 

In  this  paper,  following  the  usual  conventions,  the  x^^^  values  are  treated  as  having  constant 
values  in  the  region  of  high  transparency  between  the  bandgap  and  the  onset  of  multiphonon 
absorption,  and  the  full  x^^^  values  are  used,  which  means  the  impact  of  the  phase  matching  angle 
is  ignored.  The  requirement  to  phase  match  often  reduces  significantly  the  effective  x^^^  values 
below  those  of  the  full  values.  The  distinction  between  direct,  indirect,  or  pseudodirect  bandgaps 
is  ignored  and  the  minimum  room  temperature  gap  is  used  in  all  cases. 


A  number  of  trends  are  now  well  known  for  semiconductors.  Representative  examples  are:  a 
decreasing  thermal  conductivity  as  the  bandgap  decreases,  an  increasing  mobility  as  the  bandgap 
decreases,  an  increasing  spin-orbit  splitting  as  the  bandgap  decreases,  a  decreasing  bandgap  as  the 
lattice  constant  increases,  an  increasing  index  of  refraction  as  bandgap  decreases,  an  increasing 


dielectric  constant  as  bandgap  decreases,  and  increasing  non-linear  susceptibilities  as  the  linear  - — 

susceptibilities  increase.  The  last  three  examples  are  relevant  to  the  present  discussion  and  will  be 

□ 

described  in  some  detail  below,  followed  by  some  observations  related  to  corresponding  ab  initio  *  □ 

on 

calculations.  - — — 

The  nonlinear  optical  behavior  of  crystals  in  response  to  intense  optical  beams  is  understood  in  ?/ 
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terms  of  a  nonlinear  polarization*^’^*.  The  polarization  tensor  P(E)  is  related  to  the  electric  field 
tensor  E  through  Eq  (1)  or  Eq  (2),  where,  x®ij  =  2  dij,  K  is  the  dielectric  susceptibility,  x(2) 
(subscripts  are  omitted  subsequently)  and  are  the  nonlinear  susceptibility  coefficients. 


P(E}=K(E)*E=KoE  +  j((^)>^-\-j((3)*Ef  +  -  (1) 

Pf=>q,./*  +  ;K^2),yE/+-,  (2) 

In  Eq  (2),  the  three  dimensional  tensor  is  written  in  the  usual  "plane"  representation.  In  the 
following  discussion,  the  subscripts  will  be  omitted  except  where  helpful  for  clarity  of  meaning. 
All  non-centrosymmetric  compound  semiconductors  will  have  a  non-zero  x®  and  be  optically 
transparent.  Generally  speaking  for  infrared  applications,  all  III-V  and  II-VI  compound 
semiconductors  and  their  pseudo-ternary  analogs,  the  I-ni-Vl2  and  II-IV-V2  chalcopyrites  with 
bandgaps  of  the  order  of  2  eV  or  less,  have  large  x(2)’s,  a  number  of  compounds  from  these 
families  with  x®  's  in  the  range  of  66  pnVV  (i.e.,  AgGaSe2)  to  the  enormous  value  of  470  pmA^ 
(i.e.,  CdGeAs2)  have  received  serious  consideration  for  infrared  applications. 

The  usual  FOM  for  wavelength  conversion  efficiency  can  be  expressed  as  where  x 

is  the  effective  x-value  which  depends  on  the  specific  crystal  structure  and  the  direction  of  the 
incident  optical  beam  on  the  crystal  and  n  is  the  index  of  refraction  in  the  transparent  region.  For 
the  zinc  blende  or  the  chalcopyrite  semiconductors  the  appropriate  x  is  X^*-  If  is  traditional  to 
compare  x®  values  only  and  ignore  the  issue  of  beam  direction  in  first  order  comparisons  of  the 
nonlinear  optical  properties.  In  this  paper  trend  analysis  of  x®  vs  bandgap  and  FOM  vs  bandgap  is 
used  to  estimate  the  size  of  x®  and  the  FOM  for  wavelength  conversion  efficiency  for  known  and 
unknown  materials. 

A  trend  between  index  of  refraction  and  bandgap  was  pointed  out  by  Moss^^  in  1950.  The 


2 


Moss  expression  for  the  index  of  refraction  n,  where  the  constant  has  been  slightly  modified  by 
Ravindra  and  Srivastava^,  is  given  as  Eq  (3),  where  the  energy  in  this  and  subsequent  equations  is 
ineV: 


n  =  (108/Er  (3) 

The  data  compiled  and  modeled  by  Phillips  and  VanVechten  presented  a  corresponding 
relationship  between  the  dielectric  constant  £  or  linear  susceptibility  and  eno-gy  gap,  as  pointed 
out  by  Soref^\  as  the  Moss-like  relationships  of  Eqs  (4)  and  (5). 

e=(a/E)  (4) 

or 

/j]  1/2 

=  ihlE)  (5) 

2 

where  a  and  b  are  appropriate  constants.  Considering  that  n  =  e,  Eq  (4)  has  the  same  functional 
form  as  Eq  (3).  The  physical  basis  for  Eq  (4),  Eq  (5),  and  by  inference  Eq  (3),  is  provided  within 

29  30 

the  framework  of  the  Phillips'  bond  charge  model  ’  . 

Miller^^  in  1964  noticed  that  the  x®  values  available  for  eleven  materials  could  vary  by  several 
orders  of  magnitude  from  material  to  material.  He  found  to  his  surprise  that  the  x®  is  proportional 
to  a  product  of  the  three  linear  susceptibilities.  This  relationship  is  now  known  as  Miller's 
Empirical  Rule  and  the  proportionality  constant  is  known  as  Miller's  5ij.  Flytzanis  states  that,  in  a 
certain  sense,  this  rule  was  the  first  hint  that  nonlinear  susceptibilities  could  be  expressed  in  terms 
of  other  macroscopic  properties.  If  the  mean  value  of  Sy  determined  for  these  1 1  materials  is  used, 
then  all  of  the  d  values  can  be  calculated  to  within  a  factor  of  two.  The  rule  can  also  be  written  in 
terms  of  indices  or  dielectric  constants.  For  indices  the  rule  can  be  written  in  rationalized  mks 
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units: 


=  2(8.85 xlO^^)( (n<o-l)^)(n2a  -1)(  dij). 


(6) 


where  x®  is  in  nVV  and  na,  and  n2o)  are  the  indices  respectively  at  the  frequencies  0)  and  2®. 
Miller's  original  constant  value  for  8ij  in  his  units  was  0.38  and  the  specific  value  for  GaAs  was 
0.24.  Using  the  presently  accepted  value  of  180  pnVV  for  X36  gives  a  value  for  Miller’s  delta  of 
0.0246  in  rationalized  mks  units,  which  yields 


^2)  =0.436  (nm -if  )(n2a  -1), 


(7) 


where  the  units  of  x^^^  are  pmA^. 

Subsequently,  Flytzanis  and  Ducuing^"*,  in  a  study  of  the  specific  material-to-material  variation 
of  5ij  within  the  III-V  semiconductor  family,  found  that  Miller's  delta  correlates  remarkably  well 
with  the  dipole  moment  of  the  cation-anion  bond.  Ultimately,  using  the  bond  charge  model, 
Levine^^  provided  an  approach  for  calculating  the  value  of  6ij  for  a  given  material.  Levine's  value 
for  used  in  conjunction  with  Miller's  Empirical  Rule  overestimates  the  experimental  values  by 
20%  or  less  for  a  large  number  of  semiconductors.  This  modification  based  on  physical  insight 
greatly  improved  the  accuracy  of  Miller's  Empirical  Rule.  However,  it  was  Miller's  Empirical 
Rule  that  provided  the  insight  for  the  establishment  of  the  underlying  physical  basis  for  nonlinear 
optical  susceptibility . 

Ignoring  dispersion.  Miller's  Empirical  Rule,  as  given  by  Eq  7,  reduces  approximately,  for 
large  index  materials,  for  n^  »1,  to : 
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=(0.436)  n(8ij). 


(8) 


Combining  Miller's  rule,  Eq  (6),  and  the  Moss  relationship,  Eq  (3),  indicates  that  for 

3/2 

increases  rapidily  as  E  decreases,  going  as  (1/E)  .  In  the  same  spirit,  the  FOM  also  increases 

9/4 

rapidly  as  E  decreases,  going  as  (1/E)  .  Tabulations  of  x®  values  versus  transparency  range 
given  by  Prokhorov^^,  or  of  the  FOM  for  wavelength  conversion  versus  transparency  range  given 
by  Shay  and  Wemick  for  representative  materials,  also  suggest  the  trend  that  nonlinear  properties 
increase  rapidly  with  decreasing  bandgap.  The  prior  successes  of  empirical  property  trend  analysis 
have  served  as  the  motivation  for  the  FOM  trend  analysis  to  follow. 

The  alternate  approach  of  ab  initio  calculations  of  the  corresponding  quantities,  particularly  for 
even-order  susceptibilities  like  x®,  is  presently  not  generally  feasible.  Flytzanis^^  has  reviewed 
efforts  to  produce  approximate  theoretical  scaling  laws  for  nonlinear  susceptibilities  and  Levine  et 
al.^^  and  Zhong  et  al.'^  have  made  "near"  ab  initio  calculations  of  the  nonlinear  susceptibilities. 
These  two  approaches  are  discussed  below. 

Flytzanis  has  shown  that  odd  order  nonlinear  susceptibilities,  like  the  linear  one  can  be  cast  in 
the  form  of  scaling  laws,  essentially  power  laws  of  an  effective  parameter  which  turns  out  to  be  a 
measure  of  the  valence  electron  delocalization  and  can  be  expressed  in  terms  of  measurable 
macroscopic  quantities.  It  is  found  that  for  three-dimensional  systems  the  odd  order 
susceptibilities  scale  as: 


(2m.l) 

X 


2  (3m-l) 

=  (P  )(1/E  ), 


(9) 


where  P  is  the  average  transition-dipole-moment  matrix  element  and  m  =  0,  1,2, ....  For  x^^^ 
and  x^^^  Eq  (9)  gives 


=P^/E^ 


(m=l) 


(10) 
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and 


«1  2  5 

=(P){1IE). 


(m  =  2) 


(11) 


Note  that  the  ratio  of  /(^Ho  yields Eq  (10)  may  be  compared  to  Eq  (5),  implying 
that  p2  scales  as  E  for  scaling  as  E  and  that  scales  as  in  Eq  (12)  below: 

=(c)(llE)\  (12) 

where  c  is  an  appropriate  constant.  These  scaling  laws  were  derived  in  the  one-electron  approach, 
ignoring  electron  correlations  via  a  critical  point  analysis.  These  scaling  laws  might  become  invalid 
as  the  average  atomic  number  of  the  compound  increases  to  the  point  that  spin-orbit-splitting  is  a 
significant  portion  of  the  gap.  This  occurs  for  very  narrow  gap  materials. 

Sheik-Bahae  et  al.^^  used  a  two  band  model  to  calculate  the  scaling  of  and  found  that 
empirically  an  E  ^  dependence  is  successful  in  predicting  the  x^^^  for  a  wide  variety  of 
semiconductors  and  insulators  over  four  orders  of  magnitude.  This  is  a  somewhat  slower  but 

-6.5 

Similar  dependence  to  the  E  dependence  given  by  Eq  (12). 

33 

It  is  not  possible  to  provide  a  similar  approximate  theoretical  expression  for  even-order 
....  (2) 

susceptibilities  such  as  x  because  it  is  sensitive  to  local  properties  of  the  electron  density,  and 
integration  over  the  entire  Brillouin  Zone  is  required.  The  physical  basis  is  that  the  simultaneous 
requirement  for  large  charge  asymmetry  and  large  charge  delocalization  are  competing  processes 
resulting  in  some  cancellation  effects.  However,  the  empirical  trend  for  x^^^  and  the  approximate 
theoretical  scaling  laws  for  odd  order  susceptibilities  indicate  that,  generally,  susceptibilities  are  a 
power  law  function  of  the  bandgap,  increasing  rapidly  as  the  bandgap  decreases. 

The  "nearly"  ab  initio  calculations  of  the  nonlinear  susceptibilities  for  some  III-V  compounds 
and  Se  have  been  reported  ’  .  At  present  ab  initio  calculations  do  not  accurately  predict  the 
bandgap,  tending  to  produce  a  value  about  1.0  eV  too  small  and  a  lattice  constant  0.100  to  0.200 
nanometers  too  small.  Unfortunately,  the  value  of  x^^^  is  extremely  sensitive  to  lattice 
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constant/bond  length  and  bandgap  .  These  limitations  are  overcome  by  adding  a  constant  self¬ 
energy  correction  to  the  Hamiltonian  to  fix  the  bandgap  and  using  the  experimental  lattice 
parameters.  This  "nearly"  ab  initio  approach  then  gives  good  agreement  with  experimental  results 
for  compounds  of  elements  in  and  above  Period  IV.  However,  for  elements  and  compounds  in 
Period  V  or  greater  such  as  Te,  CdGeAs2  ,  and  InSb,  this  approach  does  not  apply,  as  once  again 
spin-orbit-splitting  has  not  yet  been  taken  into  account  The  results  for  the  five  materials  for  which 
was  calculated  also  shows  the  general  trend  that  x^^^  increases  rapidly  with  decreasing 
bandgap.  Fitting  the  data  indicates  that  x  ^  increases  as  E'  and  FOM  as  E  . 

Flytzanis^^  points  out  that  the  searches  for  scaling  laws  and  empirical  trend  analysis  for 
susceptibilities,  as  for  other  physical  properties,  serves  many  purposes,  fundamental  as  well  as 
practical.  They  allow  the  prediction  of  the  properties  of  new  materials  by  simple  inspection,  they 
introduce  economy  in  comparing  materials,  and  they  help  formulate  guidelines  for  the  search  for 
new  materials,  as  well  as  provide  guidance  for  the  development  of  theory  by  singling  out  dominant 
mechanisms.  In  the  general  case,  scaling  laws  and  empirical  trend  analysis  are  the  only  tools 
available  for  guiding  new  materials  development. 

RESULTS  AND  DISCUSSION 


Using  a  database  compiled  from  the  literature  for  a  number  of  compounds  of  interest  for 
nonlinear  optical  applications,  a  plot  of  the  second  order  nonlinear  coefficient,  x®,  versus  bandgap 
has  been  constructed,  and  is  given  in  Figure  1.  A  similar  plot  of  the  FOM  versus  bandgap  is  given 
in  Figure  2.  Where  known  or  observed,  uncertainty  for  these  values  was  noted;  otherwise  a 
minimum  of  ±20%  relative  error  was  assumed,  a  value  most  likely  too  small  an  error  for  much  of 
the  data  in  view  of  the  wide  variability  observed  from  the  literature  where  some  error  ranges 
approach  factors  of  two.  In  converting  from  relative  to  absolute  values  the  value  of  x^^^36=  180 
pmA^  from  reference  14  has  been  used  following  general  practice.  We  note  that  many  values 
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reported  in  the  literature  have  been  converted  using  much  higher  values  ranging  from  139  to  151, 
which  give  values  larger  by  50%  to  70%.  In  Table  1,  Soref  s  values  for  IV-IV  compounds  are 
used.  These  are  calculated  values  assuming  Miller’s  rule  and  using  Sy  for  SiC.  This  approach,  in 
view  of  the  saturation  shown  in  Figure  1,  may  overestimate  significantly  the  x^^^’s  of  the  narrow 
bandgap  members  of  the  family  such  as  GeSn. 

The  general  trend  of  both  and  FOM  is  to  increase  rapidly  with  decreasing  bandgap.  In 
Table  1  are  listed  the  data  (and  sources)  used  in  the  plots.  Examination  of  the  curves  shows  that  a 
single  power  law  expression  cannot  represent  the  data  well  over  the  total  bandgap  range. 
Regression  fits  were  computed  which  relate  x^^^  to  bandgap  and  FOM  to  bandgap  for  a  single  curve 
over  the  full  range  of  E,  and  two  curves.  The  least  squares  fits  are  listed  in  Table  2.  A  power  law 
was  chosen  because  it  is  consistent  with  the  results  of  the  ones  discussed  in  the  introduction, 
although  X  also  can  be  fit  reasonably  well  an  e  type  law  (x  =  1 134e'  ).  Clearly,  the  single 

curve  overestimates  large  band  gap  compounds,  while  underestimating  small  band  gap 
compounds.  For  those  materials  with  gaps  from  about  0.8  to  about  3  eV,  the  single  curve  fit 
provides  an  order  of  magnitude  estimate.  The  two  curves  provide  better  estimates  for  the  ranges  0 
to  about  1  eV  and  1  to  about  8  eV.  We  note  that  the  Miller’s  5  values  are  rather  widespread  for  this 
set  of  compounds,  covering  a  range  of  about  a  factor  of  two,  and  having  an  average  value  of  about 
0.05  m^/C. 

Calculated  values  based  on  Table  2  and  actual  values  are  listed  in  Table  3  for  several  well 
known  materials.  Table  3  indicates  that  predictions  are  within  a  factor  of  two  in  the  IR  for  both  x^^ 
and  FOM,  and  over  the  full  range  all  agree  approximately  within  the  goal  of  a  factor  of  2-5.  Further 
comparisons  are  presented  in  Tables  4  and  5.  In  Table  4,  predicted  values  and  actual  values  for 
families  are  compared  by  calculating  ratios  for  three  groupings,  using  the  lowest  value  of  the  group 
as  a  reference.  The  results  are  presented  in  columns  5  and  6,  and  they  suggest  that  the  predicted 
values  are  close  to  the  actual  values  where  known.  Again  the  accuracy  is  within  a  factor  of  two  in 
the  IR.  In  Table  5  a  similar  comparison  is  made  for  FOM.  The  results  indicate  reasonable 
agreement  except  for  the  narrow  band,  where  the  predicted  value  severely  underestimates  the  actual 
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value,  although  the  agreement  with  the  relative  value  is  within  a  factor  of  two.  Also,  we  note  that 
predicted  values  of  FOM  may  be  in  error  by  factors  more  than  two,  because  x  enters  FOM  as  the 
square.  The  trend,  however,  is  valid  for  FOM. 

The  fits  in  Table  2  are  applied  to  four  promising  relatively  unknown  materials;  1)  the  recendy 
discovered  ordered  ternary  III-V  compound  GaInP2,  predicted  by  Wei  and  Zunger^^  for  which  no 
nonlinear  properties  are  available,  2)  the  less  well  known  AgGaTe2  whose  nonlinear  properties 
have  not  been  reported  and  which  is  a  member  of  the  compound  family  whose  two  well  known 
members  are  AgGaS2  and  AgGaSe2, 3)  several  new  members  of  the  II-III2-VI4  family,  defect  or 
pseudo  chalcopyrites,  which  includes  the  well  known  mercury  thiogallate,  HgGa2S4,  and  4)  TeSe 
solid  solution  alloys  which  have  not  received  the  consideration  they  deserve  as  very  high  efficiency 
materials.  These  Equations  are  also  applied  to  three  well  known  materials,  GaAs,  ZnGeP2  and 
CdGeAs2  as  baseline  comparisons.  All  of  these  results  are  presented  in  Table  5.  As  a  final 
discussion,  several  additonal  subtrends  based  on  families  of  compounds  are  discussed. 

Random  ternary  alloy  ni-V  semiconductors  have  received  little  attention,  although  their 
can  be  very  large,  because  they  are  isotropic,  and,  therefore,  do  not  possess  birefringence  which  is 
necessary  for  conventional  wavelength  conversion  processes.  However,  in  the  ordered  state  it 
would  be  expected  that  ternaries  like  GaInP2  would  be  birefringent  uniaxial  crystals.  Stenger  et 
al."*^  have  experimentally  verified  this  expectation  of  adequate  birefringence  by  variable  angle 
spectral  ellipsometry  measurements  for  CuPt-like  ordering  in  GaInP2.  The  estimated  x®  from 
Table  4  for  ordered  GaInP2  is  130  pmA^,  a  value  large  in  practical  terms.  This  seems  to  be  a  very 
reasonable  value,  since  it  is  somewhat  less  than  the  trend  value  for  GaAs  of  214  pnVV,  as  would 
be  expected  from  its  slightly  larger  bandgap,  and  the  trend  result  for  GaAs  is  in  fairly  good 
agreement  with  the  literature*^  value  of  180  pmA^.  It  would  seem  that  films  of  ordered  GaInP2 
would  be  exceptional  waveguide  materials. 

Some  years  ago.  Bell  et  al.^^  investigated  the  properties  of  the  AgGaTe2  compounds  and 
commented  that  AgGaTe2  had  good  mechanical  properties  and  should  have  a  sufficiently  large 
birefringence  to  be  useful  for  nonlinear  optics,  although  there  had  not  been  any  studies  of  this 
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aspect  of  the  AgGaTe2  compounds.  Because  x®  data  for  AgGaTe2  is  not  available  from  the 
literature,  the  estimated  value  from  Table  4  is  233  pmA^.  Comparing  trend  results  with  literature 
values  for  AgGaS2  and  AgGaSe2  indicates  that,  for  this  family,  the  trend  result  overestimates  the 
literature  values  but  does  accurately  reflect  the  relative  values  for  these  two  materials  as  shown  in 
column  VII  of  Table  4.  In  Figure  (1),  the  II- VI  compounds  and  their  pseudo-ternaries  lie  generally 
lower  while  following  the  same  slope  which  is  consistent  with  the  actual  materials  dependence  of 
Mills’s  5.  Relying  on  the  more  reliable  relative  value  would  give  a  x®  ©f  101  pnVV  which  exceeds 
the  value  for  AgGaSe2  by  60%  and  more  importantly  a  factor  of  3.5  times  greater  FOM.  This 
material  is  certainly  worth  re-exploring. 

Radautsan  et  al.'^  have  surveyed  a  number  of  the  11-1112- VI4  compounds  and  indicate  that  their 
nonlinear  properties  are  interesting.  The  thiogallates  HgGa2S4,  HgGa2Se4,  HgIn2Se4  are 
compared  using  Table  2  and  the  bandgaps  reported  by  Radautsan  et  al.  The  trend  result  of  36.2 
pnVV  for  for  HgGa2S4  is  in  reasonable  agreement  with  the  literature  value**  for  which 
ranges  from  53.6  pm/V  to  70.4  pnVV.  So  the  trend  results  listed  in  Table  (3)  should  be  reasonable 
estimates.  For  HgIn2Se4,  the  estimate  for  x®  is  an  impressive  196  pnVV,  its  FOM  is  a  factor  of  15 
larger  than  that  of  HgGaS4!  Its  range  of  transparency  should  be  quite  large,  given  its  calculated 
range  is  from  0.86  microns  to  much  greater  than  12  microns,  the  long- wavelength  cut-off  of 
HgGa2S4. 

Tellurium  has  the  second  highest  reported  x^^^  for  a  semiconductor  and  it  is  transparent  in  the 
infrared.  However,  in  CO2  laser  doubling  experiments  it  has  been  found  that  two  photon 
absorption  and  intraband  scattering  of  free  carriers  limits  its  application  to  low  power  levels.  This 
is  a  result  of  its  narrow  bandgap.  This  problem  could  be  alleviated  by  alloying  with  Se  to  increase 
the  bandgap.  Pure  tellurium  crystals  are  quite  soft  and  deformation  induced  defects  also  limit  its 
performance'*^.  Alloy  hardening  may  adequately  improve  its  mechanical  properties.  The  optical 
properties  of  Te  crystals  alloyed  with  respectively  2%,  2.7%,  and  11%  have  been  reported^’^^. 
There  is  evidence  that  a  solid  solution  exists  for  all  concentrations'*^  where  the  lattice  constant 
decreases  in  a  linear  fashion  as  selenium  is  added  to  tellurium.  The  bandgaps  in  eV  and  the 
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respective  short  wavelength  cut-offs  in  microns  for  the  0%,  2.7%,  11%,  and  100%  cases  are 
respectively  0.358/3.46,  0.375/3.3,  approximately  0.4/3. 1,  and  1.7/0.72.  Alloying  provides  a 
method  of  continuously  varying  the  cut-off  from  3.46  to  .72  microns.  The  and  the  FOM  will 
vary  with  composition  as  the  bandgap  varies  as  indicated  by  Eq  (4)  and  (5),  but  this  dependence 
may  be  nonlinear,  as  pointed  out  by  Stuke  and  Keller"**.  Thus,  a  linear  interpolation  of  band  gap 
values  from  pure  Te  and  pure  Se  based  on  composition  will  yield  incorrect  results.  Using  the 
bandgap  composition  results  of  Stuke  and  Keller,  for  a  1  eV  bandgap/1.24  micron  cut-off  alloy, 
the  composition  Te286  is  required.  Once  a  bandgap  value  is  estimated,  x^^  FOM  can  be 
read  from  the  plots;  in  this  case,  this  compound  has  a  x^^^  of  about  700  pm/V  and  FOM  of  about 
3000pm^/vl 

The  intraband  scattering  mechanism  mentioned  above  which  limits  the  performance  of  Te  will 
be  much  less  severe  in  the  alloys.  As  the  bandgap  is  widened  the  valence  band  splitting  is 
significantly  smaller  in  comparison,  and  this  scattering  mechanism  is  much  less  important.  In 
addition,  Dubinskaya^^  reports  an  approach  which  may  overcome  this  limitation.  Te  and,  by 
analogy,  these  alloys  are  easily  bleached  to  improve  their  transparency  in  the  infrared  by  exposing 
the  crystals  to  weak  visible  light  sources.  The  impressively  large  values  for  the  nonlinear 
parameters  indicate  that  a  serious  exploration  of  this  alloy  system  is  warranted. 

The  power  law  dependencies  for  FOM  and  x^^^  ^tre  summarized  in  Table  6  and  discussed 
below.  Empirically,  the  FOM  depends  more  strongly  on  E  than  x^^^  does.  As  a  result,  reducing  the 
bandgap  is  significantly  more  beneficial  than  the  x^^^  trend  alone  might  indicate.  From  the  fit  for 

in  the  bandgap  range  >  1  eV,  and  from  the  "nearly"  ab  initio  calculation  for  the  1.4  -  2.45  eV 
range,  the  bandgap  dependences  are  similar,  as  m  =2.6  and  3,  respectively.  For  the  FOM  the 
"nearly"  ab  initio  result  and  our  results  give  respectively  m  =  6  and  m  =4.1  (using  our  large  gap 
value)  or  m  =  0.9  using  our  small  gap  value,  indicating  that  the  field  dependence  of  FOM  is 
overestimated  by  the  calculated  results. 

For  bandgaps  <1  eV,  the  bandgap  dependence  for  both  FOM  and  are  strikingly  lower  as 
saturation  sets  in.  A  comparison  of  "nearly"  ab  initio  results  for  the  narrow  bandgap  materials  has 
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not  been  made,  because  the  model  does  not  include  spin-orbit-splitting  and  it  is  useful  only  for 
compounds  made  up  of  low  atomic  number  atoms^^.  The  bond  charge  model,  the  approximate  ab 
initio  model,  and  the  "nearly"  ab  initio  model  fail  to  predict  the  saturation  observed  for  FOM  and 
although  Miller’s  rule  does,  as  shown  in  Figure  1.  However,  Miller’s  rule  does  poorly  in  the 
wide  bandgap  regime,  tending  to  overestimate  significantly.  At  the  present  time  the  bandgap 
dependence  over  the  full  range  is  treated  successfully  only  by  trend  analysis,  which  covers 
materials  with  a  gap  of  less  than  one  eV. 

A  more  careful  look  at  the  trends  in  Fig.  (1)  indicates  that  the  data  can  be  grouped  in  a  number 
of  ways  to  identify  sub-trends.  Compounds  with  Group  V  or  VI  anions  can  be  grouped  as  the 
upper  or  high  bandgap  group  while  oxides  and  complex  structures  as  the  lower  or  low  bandgap 
group.  Using  these  groupings,  there  are  at  least  two  curves  present  ~  one  for  narrow  bandgap  and 
one  for  large  bandgap.  Several  compound  families  show  similar  slopes  when  the  anion  element  is 
considered,  for  example,  {ZnS,  ZnSe,  ZnTe},  {CdS,  CdSe,  CdTe},  {GaP,  GaAs,  GaSb), 
{CuGaS2,  CuGaSe2},  {Se,  Te}  and  {AgGaS2,  AgGaSe2}  have  similar  slopes.  The  sulfides  and 
the  selenides  tend  to  fall  below  the  trend  line,  for  the  III-V  compounds,  in  agreement  with  their 
smaller  Miller’s  6’s.  The  n-IV-V2  family  members,  on  the  other  hand,  tend  to  fall  above  the  trend 
line  for  the  III-V  compounds,  in  agreement  with  their  large  Miller's  5’s.  The  {CdGeP2, 
CdGeAs2}  line  has  a  lower  slope  than  the  others,  as  does  {ZnS,  CdS,  HgS}.  The  (InP,  InAs, 
InSb}  trio  does  not  lie  on  a  line,  suggesting  that  the  data  may  be  suspect,  probably  for  the  InAs, 
because  the  line  for  InP  and  InSb  has  a  slope  similar  to  that  for  similar  groupings.  The  predicted 
value  for  In  As  is  about  1650  pm/V  ±  500.  If  the  {CdGeP2,  CdGeAs2}  data  is  correct,  this 
suggests  that  the  value  for  ZnGeAs2  should  be  about  356  pnVV  ±  75  based  on  the  direct 
bandgap  value  of  about  1.15  eV.  However,  it  is  anticipated  that  ZnGeAs2  will  not  possess  an 
adequate  birefringence,  because  its  tetragonal  distortion  is  zero  (i.e.,  c/a  =  2.0). 
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CONCLUSIONS 


As  a  means  for  quickly  estimating  the  value  of  a  material,  these  curves  offer  a  simple  and  direct 
means  for  classifying  nonlinear  optical  materials.  The  plots  show  a  strong  trend  in  the  data,  i.e., 
X®  increasing  rapidly  with  decreasing  bandgap.  The  FOM  depends  more  strongly  on  E  than 
does,  and,  as  a  result,  reducing  the  bandgap  os  significantly  more  beneficial  than  the  trend 
alone  might  indicate.  The  bandgap  dependence  of  the  “nearly”  ab  initio  and  that  of  our  empirical 
result  for  the  bandgap  near  and  above  one  eV  are  identical.  A  saturation  of  x^^  at  band  gaps  below 
one  eV  is  observed  in  our  empirical  trend,  although  the  various  models  described  fail  to  predict  this 
behavior,  and  trend  analysis  is  superior  to  Miller’s  rule  for  large  bandgaps.  This  relationship  alone 
is  of  value  in  estimating  for  a  compound  based  on  band  gap  value.  Thus,  a  non- 
centrosymmetric  compound  can  be  evaluated  in  broad  terms  as  to  its  utility  as  a  nonlinear 
optical  material.  Trend  analysis  indicates  that  films  of  ordered  GaInP2  would  be  exceptional  as  E- 
O  waveguide  materials,  that  the  FOM  of  AgGaTe2  is  an  impressive  factor  of  3.5  times  greater  than 
that  of  AgGaSe2  and  crystals  of  HgGa2Se4,  and  that  TexSe(i-x)  alloys  are  of  distinct  interest, 
because  they  are  vay  efficient  wavelength  conversion  materials  for  infrared  application. 
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Table  1.  Data  used  to  plot  against  band  gap.  Sources  of  the  data  are  given  in 
the  references  listed  at  the  bottom  of  the  table.  Sy  =  Miller’s  delta;  Index  =  index 
of  refraction.  Units  are  Gap  in  eV,  %  ***  pm/V,  FOM  in  (pm/V)2  and  5^  in  /C. 


Comoound 

Ref  No. 

Gap 

FOM 

Index 

6„ 

1 

InSb 

a,b 

0.23 

3268.0 

43322.43 

3.95 

0.056 

2 

Te 

c 

0.33 

1300.0 

3834.71 

4.794 

0.007 

3 

GeSn 

a,p,q 

0.36 

2308.0 

14614.17 

4.5 

0.02 

4 

InAs 

a,d 

0.36 

838.0 

4388.84 

3.42 

0.033 

5 

CdGeAs2 

c,e,f 

0.57 

820.0 

3821.59 

3.5 

0.029 

6 

GaSb 

d 

0.72 

1256.0 

7785.99 

3.7 

0.028 

7 

SiSn 

a,p,q 

0.84 

1010.0 

4138.01 

3.95 

0.02 

8 

SiGe 

a,p,q 

0.9 

674.0 

2224.02 

3.71 

0.02 

9 

SnC 

a,p,q 

1.2 

556.0 

1656.46 

3.6 

0.02 

10 

AgInSe2 

c,e,g,h 

1.2 

100.0 

135.87 

2.6 

0.031 

11 

InSe 

a 

1.25 

200.0 

131.19 

4.2 

- 

12 

InP 

d 

1.35 

287.0 

505.86 

3.4 

0.008 

13 

GaAs 

b,d 

1.4 

270.0 

521.22 

3.27 

0.011 

14 

CdTe 

a,c,i 

1.5 

336.0 

1551.50 

2.63 

0.077 

15 

CuInS2 

c,e,j,k 

1.53 

19.2 

5.43 

2.6 

0.007 

16 

CuGaSe2 

c,e,h,j,k 

1.7 

70.0 

60.87 

2.7 

0.020 

17 

Se 

b,c 

1.7 

159.2 

344.36 

2.64 

0.042 

18 

CdGeP2 

c,e,f 

1.72 

320.0 

788.62 

3.2 

0.025 

19 

ZnSiAs2 

c,e,f 

1.74 

180.0 

244.89 

3.2 

0.016 

20 

AgGaSe2 

c,e,g,h,k 

1.8 

95.0 

126-90 

2.6 

0.030 

21 

CdSe 

c 

1.8 

104.0 

183.87 

2.5 

0.047 

22 

Ag3SbS3 

c 

1.93 

28.0 

9.53 

2.7 

0.006 

23 

Ag3AsS3 

a 

2 

50.0 

31.75 

2.7 

0.012 

17 


24 

GaSe 

c 

2.021 

128.0 

186.59 

2.8 

0.029 

25 

ZnGeP2 

c,e,f 

2.05 

150.0 

185.21 

3.1 

0.020 

26 

HgS 

b 

2.1 

100.0 

142.24 

2.6 

0.029 

27 

GeC 

a,p,q 

2.1 

76.0 

75.02 

2.68 

0.02 

28 

AgAsS2 

a 

2.14 

50.0 

40.00 

2.5 

0.022 

29 

b-SiC 

a,p,q 

2.26 

60.0 

51.80 

2.59 

0.005 

30 

GaP 

c,d,in 

2.3 

200.0 

335.67 

3.1 

0.020 

31 

ZnTe 

a,c,i 

2.3 

184.4 

431.89 

2.7 

0.043 

32 

CuGaS2 

c,e,g. 

2.43 

22.0 

7.93 

2.5 

0.011 

h,j,k 

33 

CdS 

c 

2.485 

88.0 

181.82 

2.2 

0.071 

34 

GaS 

a 

2.5 

279.0 

988.68 

2.7 

- 

35 

map 

1 

2.6 

12.6 

11.76 

1.5 

- 

36 

AgGaS2 

c,g,h,k 

2.638 

28.0 

14.18 

2.4 

0.022 

37 

ZnSe 

c,m 

2.7 

156.8 

444.63 

2.4 

0.077 

38 

Agl 

c 

2.8 

20.0 

5.50 

2.6 

0.021 

39 

CuBr 

a 

2.91 

16.0 

8.00 

2.0 

0.033 

40 

Cul 

a 

2.95 

16.0 

5.26 

2.3 

0.017 

41 

pom 

1 

3 

6.0 

2.20 

1.6 

- 

42 

CdGa2S4 

a 

3.05 

50.0 

51.37 

2.3 

0.030 

43 

CuCl 

a 

3.17 

14.0 

7.14 

1.9 

0.043 

44 

ZnO 

a,b,c 

3.3 

3.6 

0.44 

2.0 

0.008 

45 

ZnS 

c 

3.9 

61.2 

76.96 

2.3 

0.066 

46 

La03 

b 

4 

11.1 

2.67 

2.3 

- 

47 

LiNb03 

b 

4 

10.9 

4.40 

1.9 

0.004 

48 

urea 

1 

5.9 

2.0 

0.24 

1.6 

- 

49 

SiC 

n 

6 

17.2 

3.93 

2.7 

. 
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50 

AIN 

a 

6.2 

7.7 

1.85 

2.0 

- 

51 

BBO 

1 

6.3 

1.2 

0.09 

1.6 

- 

52 

KDP 

0 

7 

0.9 

0.06 

1.5 

0.030 

53 

Si02 

b 

8.4 

0.8 

0.04 

1.5 

0.015 

All  X®  and  Miller  5  data,  except  IV-IV  compounds,  are  from  the  CRC  Handook  on  Laser  Science 
and  Technology,  Vol.  Ill,  Ed.  Martin  J.  Weber,  CRC  Press,  Boca  Raton,  FL,  1986. 

a.  Editors:  0.  Madelung  and  M.  Schulz,  Landolt-Bornstein:  Numerical  Data  and  Functional 
Relationships  in  Science  and  Technology;  Group  III:  Crystal  and  Solid  State  Physics: 
Semiconductors,  Volume  22,  subvol.  a.  Springer- Verlag,  New  York,  1987. 

b.  R.  C.  Miller  and  W.  A.  Nordland,  Phys.  Rev.  B  2, 4896  (1970). 

c.  V.G.  Dmitriev,  G.G.  Gurzadyan,  and  D.N.  Nikogosyan,  Handbook  of  Nonlinear  Optical 
Crystals,  Springer  Series  in  Optical  Sciences,  Vol.  64,  Springer- Verlag  Berlin  Heidelberg,  1991. 

d.  Lockheed  Sanders-Laser  &  Electro-optics,  Technical  Data  Sheet  (Zinc  Germanium  Phosphide), 
Nov  1,  1993. 

e.  M.  Mondal,  K.  P.  Ghatak,  and  A.  K.  Das,  Acta  Physica  Hungarica  65,  43  (1989). 

f.  H.  Hahn  and  B.  Wellmann,  Naturwissenschcften  54, 42  (1967). 

g.  H.  Y.  Ueng,  and  H.  L.  Hwang  ,  Mat.  Sci.  and  Eng.  B12,  261  (1992). 

h.  M.  Mondal,  S.  Banik,  and  K.  P.  Ghatak,  J.  Low  Temp.  Phys.  74,  423  (1989). 

i.  N.  M.  Ravindra  and  V.  K.  Srivastava,  Infrared  Phys.  19, 603  (1979). 

j.  J.  E.  Jaffe  and  A.  Zunger,  Phys.  Rev.  B  29,  1882  (1984). 

k.  J.  L.  Zyskind  and  A.  Srivastava,  J.  Crys.  Growth  81,  530  (1987). 

l.  B.  R.  Pamplin,  T.  Kiyosawa,  and  K.  Masumoto,  Prog.  Crystal.  Growth  Char  act.  1,  331 
(1979). 

m.  N.  M.  Ravindra  and  V.  K.  Srivastava,  Infrared  Phys.  19,  605  (1979). 
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n.  Elgene  R.  Nichols,  John  C.  Corbin  Jr.,  Vincent  L.  Donlan,  Technical  Report  AFAL-TR-74- 
161,  pp  61-63,  July  1974.  Available  from  NTIS. 

o.  P.  Villars  and  L.  D.  Calvert,  Pearson’s  Handbook  of  Crystallographic  Data  for  Intermetallic 
Phases',  Vols.  I,  II,  III,  American  Society  for  Metals,  Metals  Park,  Ohio  44073,  1985  (first 
printing). 

p.  R.  Soref,  Appl  Opt.  31,  4627  (1992). 

q.  R.  Soref,  J.  Appl.  Phys.  72,  626  (1992). 
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Table  2.  Summary  of  Numerical  Fit  Equations  for  (pm/V)  and  for  figure  of 
merit  (FOM)  (pmA^)2  for  E  =  0-8  eV  (full  gap  range),  0-1  eV  (narrow  gap  range), 
and  1-8  eV  (wide  gap  range). 

Full  Gap  Range:  =  291E*1*^2  correlation  coef  =  0.9 

Narrow  Gap  Range:  =  663E*-*9  correlation  coef  =  .88 

Wide  Gap  Range:  x^^^  =  510E‘2-6  correlation  coef.  =  0.78 

Full  Gap  Range:  FOM  =  666E-3-0  correlation  coef.  =  0.86 

Narrow  Gap  Range:  FOM  =  2324E'l-‘*  correlation  coef.  =  0.85 

Wide  Gap  Range:  FOM  =  1760E-4-1  correlation  coef.  =  0.76 
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Table  3.  Predicted  (pm/V)  and  FOM  (pm/V)^  values  compared  with  actual 
values  for  the  three  line  fits  for  compounds  representative  of  narrow,  mid,  and 
wide  energy  gap  values. 


£ 

Full 

Narrow 

Wide 

Actual 

Compound 

0.23 

4.91E+03 

2.45E+03 

2.26E+04 

3.27E+03 

InSb  (narrow) 

1.4 

1.52E+02 

4.91E+02 

2.14E+02 

1-80E+02 

GaAs  (mid) 

7.0 

6.88E+00 

1.17E+02 

3.37E+00 

l.OOE+00 

KDP  (wide) 

FOM: 

E 

Full 

Narrow 

Wide 

Actual 

Compound 

0.23 

5.61E+04 

l.84Et04 

7.65E+05 

4...33Et04 

InSb  (narrow) 

1.4 

2.41E+02 

1.45E+03 

4.37E+02 

2.32E+02 

GaAs  (mid) 

7.0 

1.88E+00 

1.51E+02 

5.63E-01 

L.Q.0.E-.0.1 

KDP  (wide) 
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Table  4.  Comparison  of  the  calculated  and  literature  values  of  Calculated 
values  from  equations  in  Table  2. 


Compound 

E 

feV) 

X^^’-Est. 

(oml\) 

X^^Kexp** 

fDm/Vl 

Rcl.  x'"^ 
Calculated 

Exd 

Estimated/ 
Exd  Value 

GalnPj 

1  .7 

130 

- 

GaAs 

1  .4 

21  4 

180 

- 

- 

1  .1  9 

AgGaS2 

2.73 

38.3 

2  8 

1 .00 

1 .00 

1  .37 

AgGaSOz 

1  .83 

1  07 

9  5 

2.79 

3.39 

1  .13 

AgGaTOs 

1  .356 

233 

(170)* 

6.08 

- 

- 

HgGazS^ 

2.79 

36.2 

53.6,  70.4 

1.00 

• 

0.56,0.42 

HgGa2Se4 

2.1 

75.3 

- 

2.08 

- 

- 

Hgln2Se4 

1  .45 

1  96 

- 

5.41 

- 

- 

Se 

1  .7 

1  30 

1  59 

1.00 

1 .00 

0.82 

Te.286  Se.7i4 

1 

5  1  0 

(624)* 

3.92 

- 

- 

Te 

0.33 

1780 

1300 

13.69 

8.18 

1  .37 

ZnGeP2 

2 

85.4 

1  50 

1.00 

1  .00 

0.57 

CdGeAs? 

0.57 

1090 

820 

12.76 

5.47 

1  .33 

♦♦Reference  11 

♦Calculated  using  ratio,  e.g.,  233/38.3  =  6.08;  28^6.08  =  170. 
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Table  5.  Comparison  of  calculated  and  literature  values  for  the  Figure  of  Merit. 
Calculated  values  from  equations  in  Table  2. 


Compound 

E 

FOM-Calc 

tpm/Vi2_ 

FOM-Exp 

Relative 

Calculated 

FOM 

Exp. 

Estimated/ 
Exp  Value 

GaInPz 

1  .7 

1  96 

GaAs 

1  .4 

437 

232 

- 

- 

1  .88 

AgGaS2 

2.73 

27.6 

23.4 

1 .00 

1 

1.18 

AgGaSej 

1  .83 

1  44 

1  27 

5.22 

5.43 

1.13 

AgGaTez 

1  .356 

498 

(422)* 

18.04 

- 

- 

HgGa2S4 

2.79 

25.2 

- 

1  .00 

- 

- 

HgGa2Se4 

2.1 

81.7 

- 

3.24 

- 

- 

Hgln2Se4 

1  .45 

378 

- 

15.00 

- 

- 

Se 

1  .7 

196 

51  1 

1  .00 

1 

0.38 

Te.286  Se.7t4 

1 

1760 

(4589)* 

8.98 

- 

- 

Te 

0.33 

11000 

10,048 

56.12 

19.66 

1  .09 

ZnGePz 

2 

99.9 

1  85 

1  .00 

1 

0.54 

|SI9B 

■■512Q _ 

1  .34 

"‘Calculated  using  ratio 
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Table  6.  Power-law  scaling  summary,  where  m  is  the  bandgap  power  dependence. 


_ I 

FOM 

! _ 

Bandgap  Range  teV) 

Reference 

3.0 

B 

QO 

O 

this  work 

4.1 

8  -  1 

this  work 

1.4 

1  -  0 

this  work 

^ _ 
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Figure  Captions: 


Figure  1.  Plot  of  the  second  order  nonlinear  coefficient,  versus  bandgap  for  a  number  of 
chalcopyrites  and  selected  other  types  of  compounds.  See  Table  1  for  specific  data.  The  “nearly” 
ab  initio  results  (dashed  line)  and  values  calculated  using  Miller’s  empirical  rule,  as  given  by 
equation  7,  are  plotted  for  comparison,  open  boxes:  narrow  gap;  open  circles:  wide  gap;  x:  Miller’s 
rule;  filled  diamonds:  ab  initio  calculated  values;  the  solid  lines  are  fits  for  each  range;  the  dotted 
line  is  the  fit  for  ab  initio  values. 


Figure  2.  Plot  of  the  figure  of  merit  (FOM)  versus  bandgap  for  a  number  of  chalcopyrites  and 
selected  other  types  of  compounds.  The  “nearly”  ab  initio  result  (dashed  line)  is  shown  for 
comparison.  Open  boxes:  narrow  gap;  open  circles:  wide  gap;  filled  diamonds:  ab  initio  values; 
solid  lines:  fits  for  two  ranges;  dotted  line;  fit  for  ab  initio  values. 
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